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Abstract 


Since there are only limited biogeographic studies of myxomycetes in Southeast Asia, this study was conducted to assess 
the molecular diversity of Diderma hemisphaericum, a cosmopolitan myxomycetes species that abundantly occurs in most 
moist chamber collections in tropical Southeast Asia. Specimens that contain D. hemisphaericum were gathered from 
various parts of the Philippines: (1) Los Banos, (ii) Calauan, (iii) Basud, (iv) Isarog and (v) Malilipot; and Vietnam: (i) 
Thai Nguyen (11) Ha Noi and (iii) Da Nang. Using modern molecular techniques, the common barcode nuclear marker 


18S rDNA for myxomycetes was targeted for the 50 specimen used in this study. A total of 27 ribotypes (11 from the 
Philippines and 16 from Vietnam) were obtained. The Mantel test showed that genetic differentiation among ribotypes 1s 
not correlated with geographic distances. The FST index (0.083) indicates high genetic intermixing between the 
Philippines and Vietnam. Hence, it seems that the variation across regional scale gene pool is not directly influenced by 


geographic isolation but rather by environmental selection. 
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Introduction 


Our knowledge of how certain biological organism is 
distributed worldwide is considered to be the major tenet of any 
biogeographical understanding (Pimm et al, 2014). These 
instigated many debates about patterns of global distribution 
among many macroorganisms, including pine trees (Carrillo- 
Gavilan et al., 2015), lianas (DeWalt et al., 2015), and spiders 
(Boyer et al, 2016), but only limited models have been 
proposed for microorganisms. Nonetheless, understanding the 
most appropriate biogeographical hypothesis of global 
distribution is important particularly in the light of the 
undergoing estimation of the Earth’s total biodiversity (McGill 
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et al., 2015). For example in the case of myxomycetes, we can 
assume that if all myxomycetes are widely distributed, then their 
global diversity is low. On the other hand, if myxomycetes are 
restricted in terms of their geographic distribution and their 
patterns are not by merely caused by their expected habitat 
requirements, then their global diversity should be high 
(Mitchell & Meisterfeld, 2005). So far, two opposing models for 
protist biogeography are still in debate for the biogeographic 
patterns of myxomycetes (Schnittler et al., 2017a). 

The first hypothesis is known as the Everything is 
Everywhere (EiE) model. According to this model, due to their 
relatively small spore like other 
myxomycetes (or true plasmodial slime molds) are expected to 
follow the Baas-Becking model of microbial ubiquity in nature 
(Finlay, 2002). For example, the Taimyr Peninsula, located in 
the Far North of Russia, has an extremely low temperature but 


size, microorganisms, 


Cribaria violacea, a type of myxomycetes that can be found 
usually in tropical regions, has been found in this region 
(Novozhilov et al., 1999). Another example is the recent new 
records of Redefia-Santos et al. (2017) and Bernardo et al. 
(2018). They found that what seems to be an exclusively 
temperate myxomycete is reported for the first time in the 
subtropical climate of Northern Vietnam and fragmented forest 
in tropical Philippines, respectively. This model is further 
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supported by other studies that showcased rapid rates of 
dispersal for some metacommunities of protists (see studies of 
Berdjeb et al., 2018). Most of these studies showed that protists 
are capable of long distance dispersals due to several abiotic 
factors such as wind (van der Gast, 2015), volcanic eruption 
(Bass & Boegnik, 2011) or the intervention of biological 
vectors of dispersal such as birds (Wilkinson et al., 2012), and 
even human movements (Perrigo et al., 2012). However, an 
opposing second hypothesis about these general 
biogeographical patterns for protist also has been proposed. A 
popular model that refutes the former model is the model first 
proposed by Foissner (1999) which is known as Moderate 
endemicity (ME). According to this model, even if there are 
myxomycetes species that can be cosmopolitan or globally 
distributed, portion of the ca. 1000 morphospecies (Lado, 2018) 
have their own geographically restricted distribution (Foisner, 
2006). For example, geographical barriers such as mountain 
ranges (Boenigk et al., 2018), oceans (Dagamac et al., 2017a) 
may serve as a hindrance for myxomycetes to disperse and be 
distributed like in the case of plants (/xora, Mouly & Pisivin, 
2007; Nepenthes, van der Ent et al, 2015; Cinnamomum, 
Sriramavaratharajan et al., 2016) or animals (Mimus parvulus, 
Gotanda et al, 2015; Drosophila, Craddock et al, 2016; 
Tropidurus tropiduridae, de Carvalho et al, 2007), that are 
argued to have a native restriction at a local scale (Cotterill et 
al, 2008). In the arguments published by Foissner (2008), the 
ME model contradicts EiE by proposing the following 
parameters: (i) migration rates are lower in 90% of the known 
species therefore yielding at the same time low abundance rate, 
(ii) the ability of a species to get extinct will be moderate, (iii) 
the ratio of the global species pool that are found in local scale 
setting is again moderate and (iv) at least 3096 of the 
morphospecies are considered to be restricted or endemic. In 
cases of recent myxomycete studies, it seems that this model is 
more strongly favored. For instance, there is the case of the 
genus Ceratiomyxa in which two species namely Ceratiomyxa 
morchella and Ceratiomyxa sphaerosperma (Rojas et al., 2008) 
are found only in the tropical region of Central America. 
Moreover, the patchy distributional range of the unfathomable 
myxomycete species Barbeyella minutissima (Schnittler et al., 
2000) points out the exclusivity of this myxomycete species to 
temperate Picea and/or Abies type of forests. However, these 
observations are just only recently supported by Dagamac et al. 
(2017b) using sophisticated ecological community models that 
compared the diversities and distribution of tropical 
myxomycetes. In that study, a clear bifurcation between the 
Neotropical and Asian Palaeotropical forest myxomycete 
communities was illustrated due to the number of regional 
endemic species recorded for the two forests. 
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However most of this studies concentrated at the 
morphospecies level. Recent investigations that utilized partial 
18S rDNA of some cosmopolitan myxomycete populations have 
been tested to address which of the two proposed hypothesis of 
protist biogeography myxomycetes follow. As such, the 
population genetic studies using the myxomycete models 
Badhamia melanospora (Aguilar et al., 2014) and Hemitrichia 
serpula (Dagamac et al., 2017b) seems to both agree that at a 
wider global scale, geographical barriers limit the distribution of 
myxomycete species and eventually leads to a number of cryptic 
biospecies. However, at a finer regional scale like in the case of 
Southeast Asian populations, distributions of partial 18S rDNA 
variants (ribotypes) are still unexplored. Hence, this study was 
conducted to compare the population of Diderma 
hemisphaericum, another clear cut and widely distributed 
myxomycete species across the Southeast Asian tropics. 

Diderma | hemisphaericum | (Bull) Hornemann is a 
"chocolate crinkle" looking myxomycetes species that was first 
described in Flora Danica. This species 1s reported in the 
temperate countries of Europe according to the Global 
Biodiversity Information Facility (GBIF) web portal and studies 
of Chen et al. (2013). However, recent studies in Palaeotropical 
forests in Southeast Asian countries that employed moist 
chamber cultures (see studies of Redefia-Santos et al., 2017; 
Rosing et al, 2011; Ko Ko et al, 2010), shows that D. 
hemisphaericum can just not also thrive and effortlessly found in 
a tropical and subtropical climate but also expand under 
changing climate scenarios (Almadrones-Reyes & Dagamac, 
2018). This now makes this myxomycete species a good 
candidate to check the possibility of regional range 18sRDNA 
gene flow between Philippines and Vietnam. 


Materials and Methods 


DNA extraction, amplification and sequencing 

A total of 50 samples of Diderma hemisphaericum were 
obtained from the herbarium collection of the last author. From 
these collections, two set of Southeast Asian populations namely 
from the (1) Philippines and (ii) Vietnam were selected for this 
study. Following the protocol of Dagamac et al. (2017a), 
approximately 10 mature fructifications coming from the 
herbarium collections were manually transferred into a 2 ml 
sterile safe-lock eppendorf tubes with glass beads and coarse 
sand. To maintain the aseptic condition during the transfer of 
fruiting bodies, the working space was initially washed with 
70% EtOH. To further avoid cross contamination between the 
samples, the forceps were chemically sterilized by washing it 
concurrently with 1M HCl, detergent solution and 70% EtOH. 
After transferring, the samples were then immediately 
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mechanically crushed using a microtube homogenizer (Bead 
Blaster™ 24, Benchmark Scientific, USA). Next, the DNA 
extraction using Mag-Bind Plant DNA Kit (Omega Bio- tek, 
Georgia, USA) was done following the manufacturer’s protocol. 
After, the processing of the magnetic bead set-up, the samples 
was then extracted using the robot model KingFisher™ Flex 
(ThermoScientific, Germany). Then, the amplicons from the 
extracted DNA samples has been produced following the 
modified PCR condition from Shchepin et al. 2017 and by 
making the following combination for a master mix PCR 
solution: 5.64 uL of H20, 0.50 uL of colored buffer (5x Mgo- 
Buffer), 1uL of MgCl, (Bioline, 50 mM), 0.20 uL of dNTPs, 
0.25 of forward and reverse primer (S1 and Sul9Rsp, see 
Schnittler et al. 2017b), 0.16 uL of Mango- Taq DNA 
Polymerase (5U/uL) and lastly, 2 uL of diluted DNA (10ng/uL) 
was added totalling in 10 uL in volume. Next, to check if the 
amplification process was successful, visualization of the DNA 
was performed by using 1.8% agarose gel electrophoresis for 45 
minutes at 100V. To view the generated bands on the gel 
product, UV trans- illuminator (BST- 20G- 8R) was used. All 
successfully amplified DNA samples that have produced clear 
bands were separated for the purification step using enzymes 
(Exonuclease and Alkalic Phosphatase - EAP). After subjecting 
the purified amplicons (37? Celsius, for 1 hour and 85? Celsius 
for 15 minutes) the concentration was measured using 
Nanodrop at an OD of 260nm for double stranded DNA 
(Nanodrop™ Lite Spektralphotometer Thermo Scientific™) and 
was then adjusted to get the final concentration of 10 ng/uL. 

Now to perform the cycle sequencing process, the 
following mastermix combination was prepared: 1.5 uL 
sequencing buffer, 4.69 ddH;O, 4.00 uL big Dye, 4.00 uL Half 
big dye and 0.625 uL of the S1 forward and reverse primer. 
Then, the products of the PCR process were cleaned with serial 
ethanol washing (30 ul of 96% EtOH and 100 ul of 70% 
EtOH). After overnight incubation, the samples placed inside a 
96 well plate were sent to the Zoology department of the 
University of Greifswald for sequencing purposes using the 
ABI3730 sequencer. 


Quality filtering and sequence alignment. 

The qualities of the DNA sequences were initially 
validated using the sequence reader software Chromas and were 
cross checked using the Nucleotide BLAST (Basic Local 
Alignment Search Tool) Database. Sequences that gave an E- 
value of 0.00 and had a high percentage similarity to Diderma 
sequences found in the database were then used to the 
subsequent analysis. After quality filtering of the generated 
sequences, an automatic alignment algorithm using the online 
version of MAFFT (multiple alignment using fast Fourier 
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transform) software were employed using the following 
parameters set in Dagamac et al. (2017a): E- INS- i setting with 
a 2.5 gap penalty. The new alignment file of the partial 18S 
rDNA was comprised of ca. 539 positions and was then 
manually edited by checking the sequence quality using BioEdit 
v7.0.5. Initially the both ends having the first and last ca. 20 
base pairs that provided low quality sequences were shortened. 
The final alignment were then trimmed down for only 505 base 


sites which were used for further data analyses 


Data analyses 

An analysis based on the most statistically parsimonious 
gene network for ribotypes having 95% coherence in accordance 
to the set of possible outcomes based on coalescent theory 
(Hudson, 1989) was constructed for this study. The ribotype 
networks separated by mutational steps were created using the 
software TCS v.1.21 (Clement et al., 2000) which is based on 
the calculation of the probability that pairs of ribotypes are 
similar for all other combinations of ribotypes and then connects 
only the most similar ribotypes together into a network where 
their combined probability is >95%. Therefore, the resulting 
network will remove divergent ribotypes whose true genealogy 
may be concealed by homoplastic characters (Templeton & 
Sing, 1993). 

To estimate the genetic diversity of the D. hemisphaericum 
population in (i) the Philippines and (ii) Vietnam, a RAC 
(ribotype accumation curve) were constructed using the values 
generated from the default settings in EstimateS v.9 (200 
randomizations, Colwell, 2013). Next, the correlation between 
the genetic and geographic distances was also calculated by 
constructing two independent matrices: (1) matrix consisting of 
the genetic distance from the 505 bp multiple sequence 
alignment of samples was implemented in MEGA v. 7.0 with 
the following option summary: (variation estimation method = 
Bootstrap, model/method = Maximum Composite Likelihood, 
substitutions to include = Transitions + Transversions, rates 
among sites = Gamma distributed) and (ii) the matrix which 
includes the geographical distances between the collection sites 
of any two samples being compared as calculated using an Excel 
template with the Vincenty formula (Dalgleish, 2018). The two 
matrices were then subjected to a Mantel test with 999 iterations 
in the ExtraStats function in PopTools v3.2.5 (Hood, 2010). 
Afterwards, to detect the differentiation of ribotypes between the 
Philippine and Vietnam populations, an AMOVA (analysis of 
molecular variance) and genetic differentiation calculated as Chi 
-square statistics was done using Arlequin 3.5.2.2 (Excoffier & 
Lischer, 2010) with a significance based on 10000 permutations. 
To analyze the polymorphic sites that yielded the diversity of 
the ribotypes and mean number of nucleotide differences were 
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Table 1. List of samples of Diderma hemisphaericum with their corresponding location, collector, ribotype number, and 


coordinate. 
Sample Col No, 

1 A80 

2 B123 
3 B42 

4 B04 

5 C144 
6 C143 
7 TH38 
8 C84 

9 C50 
10 TH32 
11 TH42 
12 C21 
13 TH12 
14 A23 
15 B22 
16 TH41 
17 TH31 
18 TH35 
19 CA82 
20 LB81 
21 No Label 
22 LB42 
23 BAGL23 
24 CA102 
25 LB50 
26 MAAL25 
27 LB45 
28 LB101 
29 A22 
30 A65 
31 THO1 
32 TH53 
33 C26 
34 HA38 
35 TH51 
36 TH10 
37 CA68 
38 ISAL29 
39 LB70 
40 LB24 
41 CA83 
42 ISAL21 
43 DN15 
44 DN24 
45 ISALO0O2 
46 C23 
47 HA41 
48 DN84 
49 CA61 
50 CASI 
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Collector 


JCS 
JCS 
JCS 
JCS 
JCS 
JCS 
NHAD 
JCS 
JCS 
NHAD 
NHAD 
JCS 
NHAD 
JCS 
JCS 
NHAD 
NHAD 
NHAD 
JLB 
JLB 
JLB 
JLB 
NHAD 
JLB 
JLB 
NHAD 
JLB 
JLB 
JCS 
JCS 
NHAD 
NHAD 
JCS 
NHAD 
NHAD 
NHAD 
JLB 
NHAD 
JLB 
JLB 
JLB 
NHAD 
NHAD 
NHAD 
NHAD 
JCS 
NHAD 
NHAD 
JLB 
JLB 




















Region / 
City 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Calauan 
Los Banos 
Los Banos 
Los Banos 
Basud 
Calauan 
Los Banos 
Malilipot 
Los Banos 
Los Banos 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Thai Nguyen 
Hanoi 
Thai Nguyen 
Thai Nguyen 
Calauan 
Isarog 
Los Banos 
Los Banos 
Calauan 
Isarog 
Danang 
Danang 
Isarog 
Thai Nguyen 
Hanoi 
Danang 
Calauan 
Calauan 














Country 


Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Vietnam 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Philippines 
Vietnam 
Vietnam 
Philippines 
Vietnam 
Vietnam 
Vietnam 
Philippines 
Philippines 


Ribotype 


ri 
12 
12 
r2 
r2 
r2 
r2 
r3 
r3 
r4 
r5 
r6 
r7 
r8 
19 
r10 
r10 
rll 
r12 
r12 
r12 
r12 
r12 
r12 
r13 
r14 
r15 
r16 
r17 
r18 
r18 
r18 
r18 
r18 
r18 
r18 
r19 
r19 
r20 
r21 
r22 
r23 
r24 
r24 
r25 
r26 
r27 
no amplicons 
no amplicons 
no amplicons 


Latitude N 


21°35'39" 
21°35'39" 
21°35'39" 
21°35'39" 
21°35'39" 
21°35'39" 
21°34'00" 
21°35'39" 
21°35'39" 
21°34'00" 
21°34'00" 
21°35'39" 
21°34'00" 
21°35'39" 
21°35'39" 
21°34'00" 
21°34'00" 
21°34'00" 
14°05'00" 
14°10'12" 
14°10'12" 
14°10'12" 
14°3'39" 
14°05'00" 
14°10'12" 
13°18'14" 
14°10'12" 
14°10'12" 
21°35'39" 
21°35'39" 
21°34'00" 
21°34'00" 
21°35'39" 
21°55'00" 
21°34'00" 
21°34'00" 
14°05'00" 
13°39'05" 
14°10'12" 
14°10'12" 
14°4'60" 
13°39'05" 
15°55'00" 
15°55'00" 
13°39'05" 
21°35'39" 
21°55'00" 
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Longitude E 


105?50'53" 
105?50'53" 
105?50'53" 
105?50'53" 
105?50'53" 
105?50'53" 
105?41'00" 
105?50'53" 
105?50'53" 
105?41'00" 
105?41'00" 
105?50'53" 
105?41'00" 
105?50'53" 
105?50'53" 
105?41'00" 
105?41'00" 
105?41'00" 
121?15'00" 
121?14'39" 
121?14'39" 
121?14'39" 
122257750" 
121?15'00" 
121?14'39" 
123?44'25" 
121?14'39" 
121?14'39" 
105?50'53" 
105?50'53" 
105?41'00" 
105?41'00" 
105?50'53" 
105?18'00" 
105?41'00" 
105?41'00" 
121?15'0" 
123?25'38" 
121?14'39" 
121?14'39" 
121?15'00" 
123?25'38" 
107°18'00" 
107°18'00" 
123°25'38" 
105°50'53" 
105°18'00" 
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implemented using the software DnaSP v. 5.10 (Librado & 
Rojas, 2009). 
Results 

A total of 47 out of 50 of Diderma 
hemisphaericum yielded successful amplicons (94%). The 
distribution of these 47 amplicons showed 30 came from 
Vietnam and 17 came from the Philippines. From these 47 
amplicons, 27 different ribotypes were accounted (See Table 1). 
These ribotypes were distributed wherein 11 ribotypes came 
from the Philippines population and 16 ribotypes came from 
Vietnam population yielding a ribotype density (Hd) of 0.88 
and 0.91 respectively (Table 2). The ribotype accumulation 
curve (RAC) constructed for the two populations revealed 
22.4% of the Philippines and 56.9% of Vietnam were surveyed 
as shown in Fig. 3. The Cole rarefaction curve (Fig 1) showed 
that both regions have at least the same genetic diversity in 


samples 


terms of the collected samples. 


Table 2. Summary of the polymorphic site analysis comparing 
the two countries generated by the DNAsp v.5. 1: (1) Philippines 
and (ii) Vietnam. 


PHILIPPINES VIETNAM 

Number of sequences 17 30 
Number of segregating sites 98 98 
Number of ribotypes n 16 
Ribotype diversity (Hd) 0.88 0.91 
Average number of 33.07 28.05 
differences 

0.07 0.06 


Nucleotide diversity 


The genetic differentiation of the 27 ribotypes showed a 
strong significant difference (p=0.007, a= 0.01). However, to 
check if this genetic differentiation can be explained by their 
location, the correlation analysis between genetic and 
geographic distances was performed using the Mantel test for 
the 47 sequences. The result showed that the r value (0.1260) 
was within the 9596 confidence interval (- 0.0851 to 0.1820). 
This is further corroborated when the percentage component of 
the analysis of molecular variance (AMOVA, Table 3) showed 
higher proportion “within the population" (91.678) than 
"between the population" (8.322). Moreover, the FST index 
computed in Arlequin yielded a value of 0.083. A clear 
intermixing of the ribotypes between the collections from the 
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Philippines and Vietnam was also concurred when a gene 
network was constructed (Fig. 2). 


Table 3. AMOVA (Analaysis of molecular variance) displays 
the computed FST index (0.083), with (0) indicating genetic 
material sharing and (1) as no sharing. 


Source of Sum of Variance Percentage 
Variation squares component component 
Mie 44.303 1.354 8.322 

populations 

Tm. 671.229 14.916 91.678 
populations 

TOTAL 715.532 16.270 

FST index 0.083 

Discussion 


During the past few years, there have been numerous 
studies of myxomycete biodiversity in Southeast Asia, 
especially in Vietnam (131 records, see studies of Redefia- 
Santos et al., 2018; Novozhilov et al., 2018; Tran et al., 2014) 
and the Philippines (158 records, see studies of Bernardo et al., 
2018; Macabago et al., 2017; Dagamac & dela Cruz, 2015). 
Together with Thailand (145 records, Dagamac et al. 2017b; Ko 
Ko et al., 2010; Tran et al., 2008), these two countries are 
considered to be the three most myxomycete-rich countries for 
the region. However, much of the efforts conducted for the 
myxomycetes studies in the region are concentrated on 
morphological-based ecological surveys. Comparing Vietnam 
and the Philippines, the latter country has undergone during the 
last four decades (see studies of Macabago et al., 2012; 
Dagamac et al., 2012; Macabago et al., 2010; Dogma, 1975; 
Uyenco 1973) much more of these aforementioned surveys than 
the former which has a growing myxomycete studies during 
only the last five years. This status quo brought the idea to 
conduct this pilot study for the region which explores the 
diversity of a certain myxomycetes species at a genetic level to 
address if physical geographic barriers can affect the genetic 
diversity ^ across regional populations of Diderma 
hemisphaericum. 


Ribotype diversity between Vietnam and Philippines 

In terms of genetic diversity, a total of 27 ribotypes of 
Diderma hemisphaericum have been identified. However, as 
shown in the figures la and 1b, only 22.4% of total genetic 
diversity was explored from the Philippines and 56.9% in 
Vietnam. This can be explained by the limited D. 
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Figure 1. Ribotype accumulation curve based on the Chao | estimates in (A) Vietnam and (B) 
Philippines. The dark blackline shows the Chao 1 mean and the grey highlight indicates the 5% and 
95% limits. (C) A smooth ribotype accumulation curve based on Coleman rarefaction 
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Figure 2. The constructed ribotype network using TCS. The circle size indicates the number of 
sequences represented by the ribotypes. The connecting line segments indicate mutational steps 
between alleles and the small circle in between represent the hypothetical transitional ribotypes. 
Colors of the circles indicate their collection origin: red (Vietnam) and yellow (Philippines). 
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hemisphaericum samples only coming from the available 
herbarium collections at the Institute of Botany and Landscape 
Ecology, Greifswald, Germany. In spite of this, the rarefaction 
point of the ribotype accumulation curve (Fig. 1c) showed that 
both the Philippines and Vietnam have almost similar ribotype 
diversity. In a related study of Feng & Schnittler (2017), a 
ribotype accumulation curve was also constructed for the 
diversity estimation of myxomycetes of several fruiting bodies 
collected in a deciduous forest in northeastern Germany. The 
results of their survey did not also reach an asymptotic curve 
similar with what this study obtained. With the exception of the 
species of Hemitrichia calyculata found on their study, this is 
perhaps due to the possibility of several ribotypes for a single 
morphological species. Thus, Feng & Schnittler, (2017) 
suggested that molecular assessments of diversity should be 
coupled with extensive morphological approaches to create at 
least the exact diversity estimation. Nonetheless, to achieve a 
comparable genetic diversity between the two regions, more 
samples of D. hemisphaericum coming from other herbarium 
collections and extensive field surveys for both countries is 
recommended. 


Genetic differentiation of regional population opposes global 
biogeographical pattern 

Genetic differentiation is expected when comparing 
populations coming from different countries within the region. 
However, in the case of this study, genetic differentiation is 
more prominent “within the population" than “between the 
population". This is supported by the results of the AMOVA 
and chi-squared test that showed on one hand that the genetic 
variation between (interspecific) the samples of Philippines and 
Vietnam has no significant difference and on the other hand the 
genetic variations found among specimens within a specific 
region (intraspecific) are statistically significant. This simply 
means that there is an active gene flow that might be happening 
inside the regional gene pool. The FST index showed a low 
value of 0.083, indicating high genetic intermixing between the 
Philippines and Vietnam. This can also be explained by the 
intermixing of 2 specimens having r24 variants and 7 specimens 
having r18 variants from Da Nang and Northern Vietnam 
(Hanoi, and Thai Nguyen), respectively with the r19, r23, r20, 
122, and 121 from the Philippines from the gene network (Fig 
2). Since microorganisms are well known to travel long 
distances due to their relatively small sizes, long distance 
dispersal of microorganisms for protists and fungi (see studies 
of Bonito et al., 2010; Moncalvo & Buchanan, 2008) is not that 
surprising. This could be the exact same case for, Diderma 
hemisphaericum, which was already theoretically hypothesized 
to actively disperse across long distances (Kamono et al., 2009). 
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Hence, the regional geographic barriers between Philippines and 
Vietnam could not be that strong enough to hinder the regional 
gene flow. Interestingly, this is corroborated by the Mantel test 
that showed no correlation (r value = 0.1260) between genetic 
distances and geographical distances between the populations of 
the Philippines and Vietnam. This finding now seems to support 
the EiE (Everything is everywhere) model of biogeographic 
distribution for myxomycetes that states that microorganisms 
can have almost unlimited long distance dispersal but with the 
condition that the environment will be suitable for them to live 
(DeWit & Bouvier, 2006). A similar case study of a mite- 
associated fungus, Knoxdaviesia proteae, revealed that large 
distances and geographical barriers do not isolate its population 
in the Core Cape Subregion of South Africa due to some vectors 
of dispersals like for example birds (Aylward et al, 2015). 
Interestingly, this predicament now seems to contradict several 
recent studies of myxomycetes biogeographical investigations. 
Like for an instance, a study of Dagamac et al., 2017a about 
18sRDNA variants of Hemitrichia serpula, suggested a limited 
gene flow between continents caused by geographical barriers. 
Another case study about myxomycete biogeography is the 
study of Badhamia melanospora 18sRDNA variants that 
revealed an isolation of populations between the continents of 
North America, South America and the Old World (Aguilar et 
al., 2014). But to address this, both studies look at a global 
biogeographic scale which could most likely still support the 
moderate endemicity model. Perhaps, patterns for Diderma 
hemisphaericum would be clearer if a global scale population 
genetic study is conducted. 


Implications of biogeographic distribution for conservation 
Assessing the genetic diversity of plants (Mohammadi 
& Prasanna, 2003), trees (Schaberg et al., 2008), animals (Reed 
& Frankham, 2003) and microorganisms (Li & Zhang, 2012) 
have been a budding trend for conservation genetics. This 
includes the bioresources such as food, pharmaceutical drugs, 
etc. (Wang et al., 2016; Li-Chan, 2015; Giavasis, 2014). With 
the ongoing threat of climate change, many biological organisms 
can have different responses, either a species change its genetic 
properties (genetic adaptation, Somero, 2010) or change its 
behavior, morphology and physiology (phenotypic plasticity, 
Charmantier et al., 2008) to adapt to the new climate and 
environment. However, if the species could not adapt to these 
new changes, the species population will decrease until 
extinction will occur (Thomas et al., 2004). Given these facts, 
this study contributes a new knowledge on many undermined 
but ecologically important biota such as the cosmopolitan 
tropical species Diderma hemisphaericum. It was revealed that 
the distribution of this species in Southeast Asia is widely 
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distributed in terms of morphology, but in terms of a deeper 
molecular scale, the environment plays an impact to create 
possible sympatric speciation across regional gene pool. At the 
end of the day, the diversity of the genes in a clear cut 
myxomycete species seems to be not directly influence by 
geographic isolation but most likely by environmental selection. 
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